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A polyhedron is the intersection of finitely many affine halfspaces, where an affine halfspace is a set 

R-(a,P) ={i£E": (a,x) < 0} 

for some a G R n and /3 G R (here, (a, x) = Y^j=\ a j x j denotes the standard scalar product on R n ). 
Thus, every polyhedron is the set 

P^(A,b) = {x G R" : Ax < b} 

of feasible solutions to a system Ax < b of linear inequalities for some matrix A G R mxn and some 
vector b G R m . Clearly, all sets {x G R n : Ax < 6, A'x = b'} are polyhedra as well, as the system 
A'x = b' of linear equations is equivalent the system A'x < b' , — A'x < —b' of linear inequalities. 
A bounded polyhedron is called a polytope (where bounded means that there is a bound which no 
coordinate of any point in the polyhedron exceeds in absolute value). 

Polyhedra are of great importance for Operations Research, because they are not only the sets 
of feasible solutions to Linear Programs (LP), for which we have beautiful duality results and both 
practically and theoretically efficient algorithms, but even the solution of (Mixed) Integer Linear Pro- 
gramming (MILP) problems can be reduced to linear optimization problems over polyhedra. This 
relationship to a large extent forms the backbone of the extremely successful story of (Mixed) Integer 
Linear Programming and Combinatorial Optimization over the last few decades. 

In Section [Q we review those parts of the general theory of polyhedra that are most important with 
respect to optimization questions, while in Section |2] we treat concepts that are particularly relevant 
for Integer Programming. Most of the "basic polyhedral theory" today is standard textbook knowl- 
edge. In Section [3l for some of the results we provide references to the original papers. There, we 
also give pointers to proofs of the theorems mentioned in the first two sections, where we mainly 
refer to the beautiful book by Schrijver [33 ]. There are, of course, many other excellent treatments of 
the theory of polyhedra with respect to optimization questions, e.g., in the recent survey by Conforti, 
Cornuejols, and Zambelli [7 ], in the handbook articles by Schrijver [34 ] and Burkard [4j, as well as in 
the books by Nemhauser and Wolsey [301 Grotschel, Lovasz, and Schrijver EDI . Bertsimas and Weis- 
mantel Q, Cook, Cuningham, Pulleyblank, and Schrijver [9], Wolsey [40], Korte and Vygen |[26l . 
or Barvinok [1]. The books by Ziegler [41] and Griinbaum [21] are most important sources for the 
general geometric and combinatorial theory of polyhedra, in particular of polytopes. We also refer to 
the handbook article by Gritzmann and Klee [19] as well as the one by Bayer and Lee (H. 

1. The Geometry of Polyhedra 

1.1. Some Notation. We define [p] = {1, . . . ,p} and denote by R + = {a G R : a > 0} the set 
of nonnegative real numbers. A submatrix of M G R mxn is a matrix M/ j G R /xJ for some 
7^ / C [m] and ^ J C [n] formed by the rows and columns of M indexed by the elements 
of I and J, respectively. In particular, Mj j G R is the entry in row i and column j. We write 
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M I{k = Mj )[n] and M* j<7 = M [m])J , in particular, M ijVr G IR n and M*j G IR m are the z-th row and 
the j-th column of M, respectively. The kernel of M G R mxn is ker((M)) = {x G R n : Mi = O}. 
The identity matrix Id n G IR nxn has one-entries on its main diagonal and zeroes elsewhere. 

For x G IR n and J C [n], the vector formed by the components of x indexed by elements of J is 
denoted by xj G IR J . We denote by (x, y) = Y^j=i x jVj me standard scalar product of x, y G H n . 
We consider, in the context of matrix multiplication, all vectors as column vectors, and use (...)* to 
refer to the transposed matrix or vector. We denote by O G R™ the zero vector in IR™ and by <Ej G IR n 
the standard unit vector having its i-th component equal to one, all other components being zero. 

Denoting by S(M) the set of all determinants of submatrices (formed by arbitrary subsets of rows 
and columns of equal cardinality, including the empty submatrix, whose determinant is considered to 
be one) of a matrix M G M mxn , we define 

A(M) = {§ : P) q G 5(M) U (S(M)), q^ 0} , 

and, for every finite set / V C ]R n , we set A(V) = A(M), where M is any matrix whose set 
of columns is V. Clearly, for rational matrices M G (Q mxra and (finite) sets V C (Q ra we have 
A(M),A(V)CQ. 

The encoding length of a = | G (Q with p,q € ^ relatively prime is 

(a) = 1 + fIog 2 (|p| + 1)1 + [log 2 (|g| + 1)1 . 
For a rational vector v G (Q n and a rational matrix M G Q mxn , we define 

ra m n 

( V )= n + J2(v j ) and (M)= m n + ^^(M JJ ). 
j=i i=i j=i 

Moreover, we denote by (M) max the maximum encoding length of any entry in M, as well as by 
(y} max the maximum encoding length of all components of vectors in the finite set V C (Q n . 

1.2. Basics. Most important, every polyhedron P C IR™ is convex, i.e., for all x,y G P and a G 
[0, 1], we have ax + (1 — a)y G P as well. Moreover, polyhedra are topologically closed subsets 
of W 1 . 

As the solution sets to finite systems of linear inequalities, polyhedra generalize affine subspaces, 
which are the solution sets to systems of linear equations. The criterion for Ax = b not being solvable 
via the existence of some multiplier vector A G IR m with A* A = O and (A, b) ^ generalizes to 
systems of linear inequalities in the following way (where both parts of the theorem follow easily 
from each other). 

Theorem 1 (Farkas-Lemma). For each A G IR mxn and b G IR m the following hold: 

(i) Either Ax < b is solvable or there is some A G IR+ with X b A = O and (A, b) < (but not 
both). 

(ii) Either Ax = b,x > O is solvable or there is some A G IR m with )^A > O and (A, b) < 
(but not both). 

It turns out that the "crucial solutions" to systems of linear inequalities are obtained as the unique 
solutions of certain regular linear equation systems, whose components are well-known to be express- 
ible in the following way. 
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Theorem 2 (Cramer's rule). If A £ IR nxn is regular, then, for every b E IR n , the unique solution 
x E IR n with Ax = b is given by 

det(A*. i, . . . , A* j—x, b, A* 7+1, • • • , A* n ) . , , 

Xj = ■ — -—7 — ■ — for all j E n . 

3 det(A) J 

The following estimates follows from the Leibniz formula for determinants. 
Theorem 3. There is a constant C > such that the estimates 

(a) <C-n 2 ■ (M) max for all a E A(M) 

and 

(f3)<C-n 2 - (V) max for all (3 E A{V) 
hold for allM E Q mxn and for all finite sets V C Q n . 

1.3. Polyhedral and Finitely Generated Cones. A cone is a subset if C R, n with O E if and 
ay £ K for all y E if and a E IR+. A polyhedral cone is a polyhedron that is a cone, or, equivalently, 
a polyhedron P-(A O) for some A E R mxn . 

The (convex) come /iwZ/ of a subset X C IR™ is the cone 

ccone(X) = {J2 a xX ■ X' C X, \X'\ <cc,a x >0 for all x E X'} 

(with ccone(0) = {O}) of all conic combinations of the vectors in X. A cone K C ]R, n infinitely 
generated, if there is & finite set X C R™ with if = ccone(X). Every vector in a conic hull can be 
obtained by a conic combination of few generators: 

Theorem 4 (Caratheodory's Theorem, conic version). For each X C IR™ and y E ccone(X) there is 
a linearly independent subset X' C X (in particular: \X'\ < n) with y E ccone(X'). 

The Farkas-Lemma (Part (ii) of TheoremQ]) yields a separation theorem for finitely generated cones. 

Theorem 5. If y ccone(X) for the finite set X C IR™, then there is some a E IR n with 

(a, x) < < (a, y) for all x E ccone(X) 

(i.e., ccone(X) C H^(a,0), but y H^(a,0)j. 

The following result implies that every polyhedral cone is finitely generated, which is of utmost 
importance for the theory of polyhedra. 

Theorem 6. For every matrix A E R mxn , there is a finite set X C (A(A)) n with 

P-(A,0) = ccone(X). 

The polar of a cone K C R n is the convex cone 

if = {a E R n : (a, x) < for all 2; E if} . 

The polar of a finitely generated cone ccone(X) with a finite set X C ]R" obviously is the polyhedral 
cone 

(ccone(X)) = {a E IR" : (a?, a) < for all x E X} . (1) 
From Theorem [5] we also obtain 

(P^(A, 0))° = ccone{A lji5 . . . , A m ,J (2) 

for each yl E ]R mxn . 

Moreover, from Theorem [5] one deduces (ccone(X)) 00 = ccone(X), from which one finds, by 
applying Theorem [6] as well as (two times) (Q}, the following reverse statement to Theorem [6] 
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Theorem 7. For every finite set X G R n , there is a matrix A G (A(X)) mxn w/?/j 

ccone(X) =P^(A,0). 

1.4. The Fundamental Structure of Polyhedra. The homogenization of a polyhedron P- (A, b) C 
R n (with j4 G R mxn , 6 G R m ) is the polyhedral cone 

homog(P-(A b)) = {(x,f ) 6R"©R:i2;-(6<O,?>0}. 

Applying Theorem|6]to homog(P-(^4, b)) as well as Theorem|7]to the finitely generated cone 

ccone({(x, 1) G R n © R : a: G X} U {(y, 0) G R" © R : y G F}) 

for finite sets X,F C R n , one obtains the following representation theorem for polyhedra, where 

conv(X) = { V a,x :l'cl, |X'| < oo, V Q, = 1, a x > for all x G X'} 

xeX' xex' 

denotes the convex hull of a set X C R n , and 

S + T = {s + t : s G S,t G T} 
is the Minkowski sum of S,T <Z R n . 

Theorem 8 (Weyl-Minkowski Theorem). 

(i) For every A G R mxn and b G R m f/zere are finite sets X, Y C (A(A, 6)) n wtfz 

P<(A, 6) = conv(X) + ccone(Y) . (3) 

(ii) For aZ/^nite seta I,fCR" f/j«?re are ,4 G (A(X U Y)) mxn and b G (A(X U y)) m wftfe 

conv(X) + ccone(y) = P- (A, b) . (4) 

Thus, every polyhedron can be represented by outer descriptions (intersection of finitely many 
affine halfspaces) and by inner descriptions (Minkowski sum of a polytope and a finitely generated 
cone). Clearly, algebraically an outer description may contain both linear inequalities and linear 
equations. Depending on the context, one type of description of a polyhedron can be significantly 
more convenient to deal with than the other. For instance, from inner descriptions one concludes 
readily that images of polyhedra under linear maps are polyhedra as well (see Section ITT61) . In turn, 
the fact that also preimages of polyhedra under linear maps are polyhedra is easy to prove via outer 
descriptions. From outer descriptions of polyhedra one also finds immediately that intersections of 
finitely many polyhedra are polyhedra. 

A rational polyhedron is a polyhedron for which A and b, or, equivalently, X and Y, in (O and ((U) 
can be chosen to be rational matrices, vectors, and sets, respectively. Denoting, for a rational polyhe- 
dron P C R n , by (P)mlx tne smallest number a such that there is an outer description P = P- (A, b) 
of P with a rational matrix A and a rational vector b with {(A, b)) max < a, and by (F)™°° r the small- 
est number f3 such that there is an inner description P = conv(X) + ccone(Y) of P with a rational 
finite sets X, Y C Q n with (X U Y) max < (3, we obtain the following result from Theorem[3] 

Theorem 9. There is a constant C > such that 

(PYZ? <C-n 2 - <P)E? and (P)™£ < C ■ n 2 • (P)^ 
holds for every rational polyhedron P C R". 
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In particular, whenever the rational system Ax < b (with A £ Q mxn and b £ Q m ) has any 
solution, then it also has a solution whose encoding length is bounded by a polynomial in ((A, b)}. 
This shows that the linear programming feasibility problem is contained in the complexity class NP, 
and, via Theorem [T] also in coNP. (Of course, it is well-known that this problem is even solvable in 
polynomial time). 

However, the smallest possible number of inequalities Ax < b in an outer description and the 
smallest possible cardinalities of the sets X and Y in an inner description of a rational polyhedron P 
are not bounded polynomially by each other, in general. 

The characteristic cone (or recession cone) of a polyhedron P C R n is 

char(P) = {y £ R n : x + cone{y} C P for all x £ P} . 

Theorem 10. If ^ P = P-(A,b) = conv(X) + ccone(y) C R n is a polyhedron (with A G 
R mxn , b £ R m , and X, Y C R n finite), then we have 

ccone(y) = char(P) = P-(A, O) . 

The lineality space of a polyhedron P C R n is the largest linear subspace 

lineal(P) = char(P) n (- char(P)) 

contained in char(P). 

Theorem 11. For a polyhedron / P = P-(A, b) (with A £ R mxn , b £ K rn ) we have 

lineal(P) = ker(^) . 

A non-empty polyhedron 0/PC R n with lineality space lineal (P) = {O} is called pointed. 
For most purposes, it is sufficient to consider pointed polyhedra. In fact, if a polyhedron has a non- 
trivial lineality space L ^ {O} then, in most contexts, it is sufficient to investigate instead of P its 
orthogonal projection to the orthogonal complement of L, which is a pointed polyhedron. There- 
fore, subsequently we will mainly consider pointed polyhedra. For instance, all polyhedra that are 
contained in the nonnegative orthant R™ as well as all polytopes are pointed. 

1.5. Faces of Polyhedra. For a £ R, n \ {O} and (3 £ R, we denote by 

H=(a,/3) = {x £ R n : (a,x) = 0} 

the boundary hyperplane of the affine halfspace H-(a, j3). A (proper) face of a polyhedron P C 
R n is the intersection F = P n H = (a, (3) of P with the boundary hyperplane of some halfspace 
H-(a, P) 5 P containing P. The face F is said to defined by the inequality (a,x) < (3 in this 
case. Additionally, and P itself are considered (trivial) faces of P (defined by (0,x) < —1 and 
(O, x) < 0, respectively). Non-empty faces are particularly important for optimization, because they 
are the (non-empty) sets of optimal solutions to linear optimization problems over polyhedra. 

The following result (which is a generalization of Part (i) of TheoremQ]) provides a characterization 
of those inequalities that are valid for (and thus define faces of) a non-empty polyhedron. 

Theorem 12. Let A £ R mxn , b £ R"\ a £ W 1 , and /3 G R with P = P-(A, b) / 0. The inequality 
{a, x) < f3 is valid for (all x in) P if and only if there is some A G R™ with X"A = and (A, b) < f3. 

Every face F / of a polyhedron P = P-(A, b) = conv(X) + ccone(y) (with A £ R mxn , 
b £ R m , and 1,7 C R™ finite) is a polyhedron as well with 

F = conv(X n F) + ccone(Y n char(P)) . 
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Moreover, with 

I = VqA x <b(F) = {ie[m]:FQ H=(A,*, h)} , 

we have F = {x E P : Aj^x = bi}, and, conversely, for every I C [m], the set F = {x E P : 
F C H = (A,*, is a face of P with 7 C Eq j4a .< 6 (F). 

In particular, every polyhedron has finitely many faces. Partially ordered by inclusion, they form 
the face lattice C(P) of P (including and P). For each face F / of the polyhedron P we have 
lineal(P) = lineal(P). Hence, the (non-empty) faces of pointed polyhedra are pointed polyhedra as 
well. 

The maximal elements in C(P) \ {P} are called the facets of P. A face of P is a facet if and only 
if dim(P) = dim(P) — 1 holds (where the dimension dim(P) of a polyhedron is the affine dimension 
of its affine hull aff (P)). The minimal elements in C(P) \ {0} for a pointed polyhedron P are called 
the vertices of P. The vertices of a polyhedron P are the faces of dimension zero, i.e., the faces that 
contain exactly one point (which, of course, is also called vertex). 

Theorem 13. For a point v E P in a polyhedron P = P^(A,b) (with A E R mxn and b E R m j the 
following statements are pairwise equivalent. 

(1) The point v is a vertex of P. 

(2) There are no two points x,x' E P \ {v} with v E convjx, x'}. 

(3) There is some I C [m], \I\ = n such that Aj^ is regular with v = Aj^b. 

A pointed polyhedral cone K has exactly one vertex, namely O, and therefore, one is more inter- 
ested in the minimal elements of C(K) \ {0, {O}}, which are called the extreme rays of K. The 
extreme rays of a (pointed) polyhedral cone are its faces of dimension one. The one-dimensional un- 
bounded faces of a general pointed polyhedron are called its extreme rays as well, the one dimensional 
bounded faces are the edges. The edges of a pointed polyhedron P are of the form conv{u, w} with 
two vertices v / w of P, and every extreme ray of P can be written as v + R with a vertex v of P 
and an extreme ray R of char(P) (which is pointed if P is pointed). 

An outer description P = {x E R n : A = x = b = ,A-x < b-} of a polyhedron P is irredundant if 
removing any equation or any inequalitiy from the system results in a different (larger) polyhedron, 
and turning any inequality in the system into an equation results in a different (smaller) polyhedron. 

Theorem 14. Let ^ P C R n be a polyhedron, A^ G R m(1)x ™, G ^( 2 ) E R m(2)x « 

and 6( 2 ) E R m<2) with 

P C {x E IT : A«x = ,4 (2) x < 6< 2) } (5) 
such that, for all i E [m/ 2 )], we have P % H = (A( 2 \ t +, bi). 

(i) Equality in (O holds if and only if 

• aff(P) = {iEK": A^a = and 

• /or each facet F ofP, there is ani E [m 1 - 2 )] smc/j f/ja? (A^ 2 ) $ x) < bi defines the facet F 
ofP. 

(ii) If equality holds in ©, then the right-hand- side of ([5]) is an irredundant outer description 
of P if and only if 

• A^ has full row rank, and 

• the inequalities in A^x < b^ define pairwise distinct facets of P. 

An inner description P = conv(X) + ccone(Y) of a polyhedron P is irredundant if removing any 
point from X or any vector from Y results in a different (smaller) polyhedron. 
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Theorem 15. Let ^ P C IR n &e a pointed polyhedron, and let X, Y C ]R n Z?e .s^fs with 
X O P andY O char(P). 

(i) We have P = conv(X) + ccone(Y) if and only if 

• X contains all vertices of P, and 

• Y contains a nonzero vector from each extremal ray ofchar(P). 

(ii) If P = conv(X) + ccone(Y) holds, then conv(X) + ccone(Y) is an irredundant inner 
description of P if and only if 

• X is the set of vertices of P, and 

• Y contains exactly one nonzero vector from each extremal ray of char(P). 

As for a pointed polyhedron P = conv(X) + ccone(Y) (with finite sets X, Y C IR n ), c G IR n , 
and u) = max{(c, x) : x G P} we have cj = oo if (c, y) > for some j; £ F, and u; = max{(c, x) : 
x G X} otherwise, Theorem [15] implies that a bouned linear optimization problem over a pointed 
polyhedron P attains its optimum in a vertex of P. 

1.6. Projections of Polyhedra. We mentioned above that the projection P = tt{Q) of a polyhedron 
Q C M, via a linear map tt : R d IR™ with vr(y) = Ty for a matrix T G R/ 1 ^ is a polyhedron as 
well, as, for all X, Y C IR,^, we have 

7r(conv(X) + ccone(l")) = conv(7r(X)) + ccone(7r(y)) . 

We say that a face F of Q is ^-compatible if it can be defined by an inequality (T t a, y) < /3 (valid 
for Q) for some o G R n . 

Theorem 16. For a polyhedron Q C ]R d , a linear projection it : IR^ — > R", a«<i f/ze polyhedron 
P = tt(Q), the map defined via F h-> tt(F) is an isomorphism between the sublattice of C(Q) formed 
by the n -compatible faces and C{P). 

An extended formulation for a polyhedron P C ]R n is an outer description Q = P-(^4, b) of some 
polyhedron Q C ]R d along with a linear projection ir : H d — > IR n with 7r(Q) = P. Such an extended 
formulation of a polyhedron P can be much simpler (e.g., in terms of the number of inequalities) than 
any outer description of P if the complexity of the facets of P is hidden in lower dimensional parts of 
the face lattice C{Q) (see Theorem [T6l). As, due to 

max{(c, x) : x G P} = max{(T t c, y) : y G Q} (for all c G IR n ) , 

linear optimization problems over P = 7r((5) can be solved by solving linear optimization problems 
over Q, extended formulations play an important role in modern mathematical optimization (see |6l). 

We conclude the first part of the article by considering the question how to derive an outer de- 
scription of a polyhedron from an extended formulation. Here, the fundamental result follows via 
Theorem [121 

Theorem 17. Let Q = P-(D,g) C K d be a polyhedron with D G IR 9Xd and g G IR 9 , suppose the 
linear projection tt : H d — > IR n is defined via 7r(y) = Ty for a matrix T G IR nxd and all y G IR^, 
and let T be any matrix whose rows form a basis o/ker(T). 

If L G IR™ xrf is a matrix whose rows generate the projection cone 

{A G : X^DT*) = O} = ccone{L M , . . . ,L m ^} , 
then every A G K mxn with AT = LD satisfies 

7T(Q) =P^(A,b)ri7T(U d ) 

with b = Lg. 
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If, in the situation of Theorem [T7J the projection tt : 1R, d — > R n is the orthogonal projection to the 
first n coordinates, then the projection cone is simply 

{A G R q + : (D+j, A) = for all j G [d] \ [n]} , (6) 

and A can be chosen to consist of the first n columns of LD. If furthermore n = d — 1 holds, then a 
finite generating system for the projection cone © can be obtained from the relation 

r+s 

{A G K^~ s : ^2 = 0} = ccone{<ei - <B k : i G [r], k G [r + s] \ [r]} 

i=i 

for all numbers 5\, . . . , 5 r > > S r+ i, . . . , 5 r+s (which one can easily establish by induction on r + 
s). 

Theorem 18 (Fourier-Motzkin elimination). For a polyhedron Q = P-(D, g) C R d (with B G W xd 
and g G H q ) and the sets 

J>° = G [ g ] : £) i)d > 0} 

J=° = {i G [q] : A,d = 0} 
J<° = {i G [g] : £). )d < 0} , 

polyhedron in ]R, a!_1 arising from Q by orthogonal projection to the first d — 1 coordinates is the 
set of all x G IR, rf_1 f/zaf satisfy the following system: 

(A,*, a) < ft /or a// i G I =0 

(D k)d D e>ir - D i)d D kiin x) < D k)d g t - D l4 g k for all k£l >0 ,££ I <0 

By applying the Fourier-Motzkin elimination method (Theorem IT8T > iteratedly, one can compute 
from an outer descriptions of a polyhedron an outer description of its orthogonal projection to any 
coordinate subspace. Note that the sizes of the descriptions may blow up exponentially (even if, after 
each iteration, one removes redundant constraints). 

The case of general projections can be solved by Fourier-Motzin elimination as follows, where we 
adopt the notation from above. Let h be the rank of T, and assume that the submatrix X| A i r^i of T is 
regular. Denoting by T G H dxd the regular matrix that equals T^ * on its first h rows and that has the 
vectors e», for i = n + 1, . . . , d, in its other rows, and by tt : R d — > R n the orthogonal projection 
to the first h coordinates, we have Tr(y) [fl] = n{fy) for all y G R d . Thus, with D = DT' 1 and 
Q = p< (D, g), we find P = {x £ vr(]R d ) : Xfl G n(Q)}. 

One reason for the interest in computing projections of polyhedra is that one can convert inner 
to outer and outer to inner description by means of such computations. Indeed, it follows from the 
discussion of homogenization and the concept of polarity that all such conversions can be done by 
any method that computes, for a given finite set X C ]R n , a matrix A G R mxn with ccone(X) = 
P-(^4, 0) (see Theorem [7]). Denoting by T G H nxd a matrix whose set of columns is X, we have, 
by definition, ccone(X) = tt(Q), where tt : H d — > IR™ is the projection defined via Tr(y) = Ty for 
all y G H d and Q = R^. Thus, computing an outer description of tt(Q) from the outer description 
= p<( - ld d , O) of Q solves the problem. 
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2. POLYHEDRA AND INTEGRALITY 

Integral points in polyhedra play a crucial role in integer linear programming and combinatorial 
optimization. The fundamental concept here is the integer hull 

Pi = conv(P n r E n ) 

of a polyhedron P C R™. In this second part, we describe the most fundamental results on integer 
hulls, which are important, since, on the one hand, they satisfy 

max{(c, x) : x £ P H 1i n } = max{(c, x) : x G Pi} 

for all c G R", while on the other hand, as convex objects, they are much more convenient to deal with 
than the discrete sets P n 7L n . In particular we will treat integral polyhedra, i.e., polyhedra P with 
Pi = P. Identifying a polyhedron P as integral is particularly pleasant, as in this case integer linear 
optimization problems over P can be treated as (continuous) linear optimization problems over P. 

Theorem 19. For A G Q mxn and b G Q m , defining an integral polyhedron P-(A, b), one can, for 
every c G (Q n , find x* G P-(A, b) H 7L n with 

(c, x*) = max{(c, x) : x G P~(A, b) n Z n } 

(or conclude that no such x* exists) in time bounded by a polynomial in (A) + (b) + (c). 

2.1. Integer Points in Polyhedra. Similarly to the parameterization ([3]) of points in a polyhedron 
by an inner description, for a rational polyhedron P, we can also parameterize the set of all integral 
points in P. For a finite set Y C IR n , we denote by 

mono(Y) = a y y : a y G N for all y G Y} 
(with N = {0, 1, 2, . . . }) the submonoid of R n generated by Y. 

Theorem 20. For every rational polyhedron P C IR™ there are finite sets I,ycP w/f/i 
P n Z n = X + mono(Y) W char(P) = ccone(Y) , 

such that {X U Y) 

max is bounded by a polynomial in n and (P)max r ( an d thus, by a polynomial in n 

and (P)™£). 

A finite set H C % n is called an integral Hilbert basis (of ccone(-ff)) if 

ccone(-ff) Pi Z n = mono(P') 
holds. The existence statement of the following theorem is a special case of Theorem |20l 

Theorem 21. Every rational polyhedral cone K has an integral Hilbert basis; if K is pointed, then 
the integral Hilbert basis of K is uniquely determined. 

Moreover, Theorem [20] implies that, whenever the rational system Ax < b (with A G Q mxn 
and b G Q m ) has any integral solution, then it also has an integral solution whose encoding length 
is bounded by a polynomial in ((A, b)). This shows that the integer linear programming feasibility 
problem is contained in the complexity class NP (it is, however, not in coNP, unless NP = coNP). 

Finally, from Theorem [20] one can derive that the integral hulls of rational polyhedra are rational 
polyhedra as well. 

Theorem 22. For each rational polyhedron P C IR™, the integer hull Pi ofP is a rational polyhedron, 
for which {Pi)m&x i J bounded by a polynomial in n and (P)™^- 
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The number of facets of P\ is, however, in general not bounded polynomially in the number of facets 
of P. Furthermore, the integer hull of a non-rational polyhedron needs not even be a polyhedron (see, 
e.g., P = {(x 1 ,x 2 ) G R 2 : x 2 < V2xi,x t > 1}). 

A generalization of Theorem |22] to mixed-integer hulls of rational polyhedra holds as well: If 
P C R, n is a rational polyhedron, then conv{x G P : xj G Z J } is a rational polyhedron for every 
J C [n]. 

We end this section by a crucial criterion for integrality of (pointed rational) polyhedra. 

Theorem 23. A pointed rational polyhedron is integral (i,e„ Pi = P) if and only if all vertices of P 
are integral. 

2.2. Total Dual Integrality. A rational system Ax < b with A G Q mxn and b G Q m defining a 
non-empty polyhedron P-(^4, b) ^ is called totally dual integral (TDI) if, for every c£P with 
max{(c, x) : Ax < b} < oo, there is some y k G N m with ^y* = c and 

(6,^)=mm{(6,y):^y = c,yGR™}, 

i.e., the dual problem to max{(c, x) : Ax < 6} has an integral optimal solution y*. It is convenient to 
consider TDI all rational systems Ax < b with P-(A,b) = as well. 

Theorem 24. A system Ax < b with integral matrix A G Z mxn and rational right-hand-side vector 
b G (Q m is TDI if and only if, for every face F^0ofP- ( A, b), the set 

{A^ : i G Eq^^OF)} 

is an integral Hilbert basis. 

Applying Theorem [21] to the rational normal cones 

{c G R n : (c, x*) = max{(c, x) : x G P} for all x* G F} 

of the (minimal) faces F of a rational polyhedron P, one can construct TDI-systems as in the following 
result. 

Theorem 25. Let P C R n fte a rational polyhedron. 

(1) There is an integral matrix A G Z mxn a«<i a rational vector b G Q m w/?/j P 
such that Ax < b is TDI, where b can be chosen to be integral if P is integral. 

(2) If there is a rational matrix A G Q mxn and an integral vector b G Z m with P 
such that Ax < b is TDI, then P is integral. 

Thus, every rational polyhedron admits an outer descriptions by a TDI-system with integral (left- 
hand-side) coefficient matrix, though, in general, this will not be an irredundant description in the 
sense of Section [T31 Nevertheless, such a description can provide much structural insight, e.g., within 
the theory of cutting planes for non-integral polyhedra. The algorithmic problem, however, to decide 
whether a given rational system Ax < b is TDI is coNP-complete (even for A restricted to the set of 
node-edge incidence matrices, see Section 1231 of undirected graphs). 

Outer descriptions of (integral) polyhedra with integral right-hand-side vectors are very important, 
because due to Theorem[25ll2]) (and the definition of TDI) they yield strong duality relations for certain 
integer linear optimization problems. 

Theorem 26. If the system Ax < b is TDI with A G Q mxn and an integral right-hand side b G 7L m , 
and P = P- (A, b) ^ holds, then, for every c G 7L n with max{(c, x) : x G P} < oo, we have 

max{(c, x) : Ax < b, x G 7L n } = mm{(b, y) : A l y = c,y > Q,y G Z m } . (7) 



= ?-(A,b) 
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2.3. Total Unimodularity. A matrix A G { — 1,0, l} mxn is totally unimodular (TU) if all square 
submatrices of A have determinant —1, 0, or 1. The (version for pointed polyhedra of the) following 
result is a consequence of Theorem [23] and Theorem [13] (see part (f3]>) obtained by Cramer's rule 
(Theorem [2]). 

Theorem 27. If A G {—1,0, l} mxn is totally unimodular and b G W 1 is integral, then P-(A, b) is 
an integral polyhedron. 

If A G {—1,0, l} mx " is totally unimodular, then so are all submatrices of A, (A,ld m ), (A, —A), 
and A t , where the latter allows to conclude the following from Theorem [27] 

Theorem 28. If A € {—1, 0, l} mx ™ is totally unimodular, then, for every b G Q m , the system A < b 
is TDI. 

In particular, if A G {—1, 0, l} mxn is totally unimodular, then the strong duality relation $7} holds 
for all integral b G 7L m and c G TL n (for which the respective optimal values are finite). 

The following criterion is extremely useful for establishing total unimodularity of matrices. 

Theorem 29 (Criterion of Ghouila-Houri). A matrix A G {—1, 0, l} mxn is totally unimodular if and 
only if, for every subset I C [m] of row indices, there is a partitioning I = (with I + f)I~ = 0) 

such that 

^^,*-^A,*e{-i,o,i}" 

holds. (Clearly, a similar characterization via all subsets of column indices holds.) 

From Theorem [29] one readily deduces that every matrix with entries from { — 1, 0, 1} that has at 
most one positive and at most one negative entry per column is totally unimodular. In particular, 
the node-arc incidence matrix inc(-D) G {—1,0, 1}^ X ^ of a directed graph D = (V, A) (with A C 
V x V), defined via 



mc(D) VjCl 



— 1 if a = (v, w) for some w G V 
1 if a = (u, v) for some u G V 
otherwise 



for all v G V and a G A, is totally unimodular. Thus, the set 

{i£R A : mc(D)x = OJ<x<u} 

of circulations in D respecting the integral lower and upper bounds £, u G 7L A is an integral polytope. 

The node-edge incidence matrix inc(G) G {0,l} VxE of an undirected graph G = (V,E) is 
definied via 



inc(G% 



if e = {v, w} for some w G V 
otherwise 



for all v G V and e G E. Theorem [29] yields that the node-edge incidence matrix of a graph G = 
(V, E) is totally unimodular if and only if the graph is bipartite, i.e., there is a partitioning V = S l+J T 
(with S n T = 0) such that E C {{s, t] : s G S,t G T}. In particular, the matching polytope of a 
bipartite graph G = (V, E) (i.e., the convex hull of the characteristic vectors x{M) G {0, 1} E — with 
x(M) e = 1 if and only if e G M — of all matchings M C E — with e n e' = for all e, e' G M, 
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e ^ e') has 



x e < 1 for all u G V 



j; e > for all e € E 1 



as an outer description. 

The most important class of totally unimodular matrices is formed by the network matrices. Such a 
network matrix N G {—1,0, 1}^ tX/1 arises from a directed graph D = (V, A) and a subset C A 
of |V| — 1 arcs that, viewed as undirected edges, form a spanning tree on V, by setting, for each 

a' £ At and a = (v, w) 6 A, 



iVa',a 



1 if the w-w-path in At uses a' in forward direction 
— 1 if the v-w-path in At uses a' in backward direction . 
otherwise 



Theorem 30. Network matrices are totally unimodular. 

In fact, network matrices are the crucial building blocks of the whole class of totally unimodular 
matrices, as every totally unimodular matrix arises from networks matrices and two special totally 
unimodular 5 x 5-matrices by certain operations that preserve total unimodularity. From such struc- 
tural results one can derive polynomial time algorithms for testing matrices A for total unimodularity 
(while, as remarked at the end of Section 12.21 testing a system Ax < b for total dual integrality is 
coNP -complete) . 



3. Pointers to Literature 

Whenever possible, we provide, for the results mentioned in the article, pointers to proofs in l33l . 
which is also a great source for historical notes and many more references than listed here. 

The Farkas-Lemma (Theorem [TJ dates back to work by Farkas as well as by Minkowski at the end 
of the 19th century (see [33 , Sect. 7.3]). There are several possibilities to prove the theorem. In fact, 
the Farkas-Lemma is a special case of more general separation theorems in convex analysis (see, e.g., 
[31, Chap. 2]). A particularly nice and purely linear algebraic proof is due to Conforti, Di Summa, and 
Zambelli [8], who derive the Farkas-Lemma from the corresponding obstruction of the solvability of 
linear equation systems. Theorem[2]is due to Cramer (1750) (see fl33l Sect. 3.1]). Theorem [3] follows 
from ll33l Thm. 3.2] (see also |[T2"1 ). Theorem[4]is due to Caratheodory [5] from 1911 (see also ll33l 
Thm. 7.1]). Theorems [6l El and [8] have their origins in the work of Farkas, Minkowski l28l . and 
Weyl |39l (see also [ 33, Sect. 7.2]). The statements on the components of the vectors and on the entries 
of the matrices (as well as Theorem© follow, e.g., from (the proofs in) [33, Sect. 10.2]. An elementary 
proof of Theorem[6]can be found in [24] (see also IT331 Cor. 7.1a]). For proofs of Theorems ITOl and [TT1 
(on the characteristic cones and lineality spaces of polyhedra), we refert to lf3~3l Sect. 8.2]. Theorem[T2l 
is Cor. 7.1h in [33]- For the other statements in Section [T31 see 1331 Sect. 8.3-8.9]. Theorem [T6l is 
folklore (we are not aware of any other explicit reference, thus we refer to l25l ). Theorem [T71 usually 
is formulated for orthogonal projections to coordinate subspaces only (see, e.g., Q Sect. 2.4]). An 
explicit proof in the general setting can be found in |[25l . The Fourier-Motzkin method (Theorem IT~8l) 
is treated in [33, Sect. 12.2]. The method is due to Motzkin [29], where the idea goes back to work 
of Fourier in the early 19th century. For the algorithmic problem of converting representations of 
polyhedra, we refer to the survey by Seidel [35] and to the software system polymake by Gawrilow 
and Joswig lPT5l flhttp : / /www . opt . tu-darmstadt . de /polymake /| >. 
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The proof of Theorem [T9]relies on both that (continuous) linear programs (Khachiyan 11221 ) as well 
as systems of linear Diophantine equations can be solved in polynomial time (see IT331 Thm. 16.2]). 
For proofs of Theorems [201 |2T1 and[22l we refer to [33 ( Sect. 16.2-16.4, 17.2]. The fact that the integer 
hull of a rational polyhedron is a rational polyedron (Theorem 1221 is due to Meyer 112711 . the notion 
of Hilbert bases has been introduced by Giles and Pulleyblank ffTTl . where the ideas of the proof 
of Theorem [21] date back to Gordan lfl8ll . The concept of total dual integrality has been invented 
by Edmonds and Giles lfl3l[T4l . See IT331 Sect. 22.3] for proofs of Theorems l24l and l25l (the results 
being due to lfT7l and Schrijver ll32l ). The coNP-hardness of the TDI-property has been established 
by Ding, Feng, and Zang ifTTTl . Proofs of the results on total unimodularity can be found in ll33l 
Chap. 19]. The connection between totaly unimodular matrices and integral polyhedra (Theorem l27l 
and a similar characterization of total unimodularity) is due to Hoffman and Kruskal [23]. Theorem|29l 
has been proved by Ghouila-Houri [16]. The total unimodularity of network matrices (Theorem l30l) 
is due to Tutte |[38l . The decomposition theorem for totally unimodular matrices mentioned after 
Theorem [30] has been proved by Seymour |[36l . Cunningham and Edmonds [ 10] derived a polynomial 
time test for total unimodularity from that theorem, the asymptotically fastest known algorithm is due 
to Truemper lt37l . 
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